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ABSTRACT 

The dynamics of the Magellanic Stream (MS) as a series of clouds extending from the Magel- 
lanic Clouds (MCs) to the south Galactic pole is affected by the distribution and the amount of 
matter in the Milky Way. We calculate the gravitational effect of the Galactic disk on the MS in 
the framework of modified Newtonian dynamics(MOND) and compare with observations of the 
Stream's radial velocity. We consider the tidal force of the Galaxy, which strips material from the 
MCs to form the MS, and, using a no- halo model of the Galaxy, we ignore the effect of the drag of 
the Galactic halo on the MS. We also compare the MONDian dynamics with that in logarithmic 
and power-law dark halo models and show that the MOND theory seems plausible for describing 
the dynamics of satellite galaxies such as the MCs. Finally, we perform a maximum likelihood 
analysis to obtain the best MOND parameters for the Galactic disk. 

Subject headings: dark matter -Galaxy: halo -gravitation - Magellanic stream - MOND. 



1. Introduction 

It has long been known that Newtonian grav- 
ity is unable to describe the dynamics of galaxies 
and clusters of galaxies correctly and that there 
is a discrepancy between the visible and the dy- 
namical masses of galaxies. Nowadays, most as- 
tronomers believe that the universe is dominated 
by dark matter and that we can observe the grav- 
itational effects of this matter in large-scale struc- 
ture. Many observations have been dedicated to 
measuring the amount and nature of the dark 
matter. Gravitational microlensing is one of these 
experiments, designed to indirectly detect the 
dark matter of the halo in the form o f MACHOf] 
(|Alcock et al. 200(1 IPaczvhski 19861 ). After al- 
most two decades of microlensing experiments, it 
seems that MACHOs contribute a small fraction 
of the halo mass and that the microlensing events 



toward the Large and Small Magellanic Clouds 
are mos t probably due to self-lensing by the 
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Clouds (ISahu 20031 : iTisserand fc Milsztain 2004 ; 



Other observations, such as of the abundances 
of the light elements from nucleosynthesis in the 
early universe and recent data from the Wilkinson 
Microwave Anisotropy Probe ( WMAP) data, also 
rule out baryons as the dark matter. The best- fit 
baryon abundance, based on WMAP data, to a 
ACDM model is tl b h 2 = 0.0237t£;^, implying 
a baryon-to-photon ratio of r\ — 6.5^*3 x 10~ 5 
(jSnergel et al. 200.1 ). For thi s abundance, stan - 
dard big bang nucleosynthesis (jBurles et al. 20011 ) 
implies a primordial deuterium abundance rela- 
tive to hydrogen of [D/H] = 2.37±g;^ x lO" 5 , 
which is in the same range as observed in Lya 
clouds (Kirkman et al. 2003): analysis of QSO HS 
1243+3047 yields a D/H ratio of 2.42 
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WIMPsS, as the other candidate for the so- 
called the cold dark matter (CDM), are moti- 
vated from non-standard particle physics. The 
problem with WIMPs is their small cross section 



Weakly interactive massive particles 
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with ordinary matter, which would makes them 
difficult to observe, to date, no significant signal 
from experiments aimed at detecting such parti- 
cles has been reported dPretzl 20021: iGreen 2 004: 



Goodman fc Witten 1985t Jungman et al. 1996 ; 
Mirabolfathi 20051 ). 

Another growing approach to explain the miss- 
ing matter in the universe is a possible alternative 
theory of gravity. In these models, modification 
of the gravity law compensates for the lack of 
matter. One of the most famous alternative the- 
ories is modified newtonian dynamics (MOND), 
which was been introduced by Milgrom (1983). 
According to this phenomcnological theory, the 
flat rotation curves of spiral galaxies at large ra- 
dial distances can be explained with a modifica- 
tion of Newton's second law for accelerations be- 
low a character isti c scale of an — 1 x 10~ 8 cms~ 2 
(iMilgrom 19831: iBekenstein and Milgrom 1984 : 



Sanders fc McGaugh 20021 ). This theory, in ad- 
dition to the acceleration parameter, an, employs 
an interpolating function to connect and MON- 
Dian regimes. It has been shown that this 
simple idea may explain the dynamical motion 
of galaxies without resorting to dark matter 
dPointecouteau fc Silk 20051: iBaumgardt et al. 2005 
IScarpa et al. 2004 ). The parameters of this model 
can be fixed by various observations, such as 
the rotation curves of spiral galaxies. Recently, 
Famaey & Binney (2005) studied MOND in the 
context of the rotation curves of spiral galaxies, 
using different choices for the MOND interpolat- 
ing function. Bckenstein (2004) presented a rela- 
tivistic version of MOND, considering a Lorentz- 
covariant theory of gravity. Zhao & Famaey (2006) 
also examined different interpolating functions by 
fitting to a benchmark rotation curve, and they 
proposed a new set of interpolating functions that 
satisfy both weak and strong gravity. One of 
the main tests of MOND is its ability to produce 
power spectra that are compatible with the struc- 
tures that we observe. Sanders (2001) showed 
that perturbations on small comoving scales enter 
the MOND regime earlier than those on larger 
scales and, therefore, evolve to large overdensities 
sooner. Taking the initial power spectrum from 
the cosmic microwave background, Sanders found 
that the evolved power spectrum resembled that 
of the standard CDM universe. However, studies 
by Nusser & Pointecouteau (2006) show that for a 



spherically symmetric perturbation, in the frame- 
work of MOND, the gas temperature and density 
profiles evolve to a universal form that is indepen- 
dent of both the slope and amplitude of the initial 
density profile. While the obtained density profile 
is compatible with results from XMM — Newton 
and Chandra observations of clusters, the tem- 
perature profiles do not reconcile with the MOND 
prediction. 

Another approach to examining MOND could 
be through the dynamics of the satellite galax- 
ies around the Milky Way (MW). Read & Moore 
(2005) compared the precession of the orbit of the 
Sagittarius dwarf in a thin axisymmetric disk po- 
tential under MOND with that resulting from a 
thin disk plus nearly-spherical dark halo. They 
find that improved data on the leading arm of the 
Sagittarius dwarf, which samples the Galactic po- 
tential at large radii, could rule out MOND if the 
precession of the orbital pole can be determined to 
accuracy of a degree. Here we extend this work, 
using the dynamics of the Magellanic Stream (MS) 
to test MOND. The MS is a narrow band of neu- 
tral hydrogen clouds that extends along a curved 
path, starting from the Magellanic Clouds and ori- 
ented toward the south Galactic pole (van Kuilen- 
burg 1972; Wannier & Wrixon 1972; Mathewson 
et al. 1974). We use the MONDian gravitational 
effect of a thin Kuzmin disk to represent the lumi- 
nous part of the MW, and we obtain the dynam- 
ics of the MS, treating it as a quasi-steady flow. 
Using a maximum likelihood analysis, we find the 
best model parameters from an analysis of the MS 
radial v elocity data report ed by the Parkes Obser- 
vatory (jBriins et al. 20051 ). 

The paper is organized as follows: In fJ2J we 
give a brief review of the dark halo models for the 
Galaxy and of MOND theory as an alternative to 
dark matter. In we introduce the Magellanic 
Stream and the main theories describing the dis- 
tribution of gas along it. In 21 we compare the ob- 
served dynamics of the MS with predictions from 
CDM models and MOND theory, considering the 
MS as a quasi-steady flow. We conclude in Sj5j 

2. Dynamics of The Milky Way: MOND 
Versus CDM 

As pointed out in the previous section, there are 
two approaches for dealing with the problem of the 
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rotation curves of spiral galaxies. The first is to 
consider a dark halo around the galaxy that is 10 
times more massive than its luminous component, 
and the second is to invoke an alternative theory 
such as MOND. In the first part of this section, we 
introduce MOND theory and consider a Kuzmin 
disk for the mass distribution of the Galaxy, and in 
the second part various dark matter models for the 
halo mass distribution in the MW are presented. 

2.1. MOND 

As mentioned in §1, the shortage of luminous 
matter in spiral galaxies and large-scale structures 
has motivated some astrophysicists to apply alter- 
native theories. One possible approach to solve 
the dark matter problem is a change in Newtonian 
dynamics at small accelerations, so-called modi- 
fied Newtonian dynamics (Milgrom 1983). In this 
model, Newton's second law is altered to read 



H{- — l -)a = -V$at, 
a 



(1) 



where $jy is the Newtonian gravitational poten- 
tial, a is the acceleration, clq is an acceleration 
scale below which the dynamics deviates from the 
standard Newtonian form, and fi(x) is an interpo- 
lating function that runs smoothly from fi(x) = x 
at x -C 1 to n(x) = 1 at i > 1 (i.e., the latter 
condition recovers Newtonian dynamics at acceler- 
ations that are large compared with ao). Any suc- 
cessful underlying theory for MOND must predict 
an interpolating function. The functional form of 
fj,(x) and the value of free parameters such as ao 
are usually chosen phenomcnologically. The stan- 
dard interpolating function used to fit with obser- 
vations was introduced by Bekenstein & Milgrom 
(1984): 



(2) 



Other func tions have been suggested for n(x), such 
as those of ( Famaev fc Binnev 20051 : Bekenstein 2004 ; 



Zhao et al. 20061) . The typical value for the pa- 
rameter ao, obtained from analysis of a sample of 
spiral galaxies with high-quality rotation curves 
is ao = 1.2 ± 0.27 x 10 -8 cms~ 2 (Begeman et al. 
1991). Taking the divergence of both sides of 
equation (fTJ and substituting Poisson's equation 
in the right-hand side of equation (fT]) yields 



V • /x - a 

a 



-AnGp. 



(3) 



For a given distribution of matter on the right- 
hand side of equation ([3]), the dynamics of the 
structure, depending on the interpolating func- 
tion, can be found as 

M (iL)a = G / G(x-x>(x')dx' 3 +Vxh(x) (4) 
ao J 

where G(x) is the Green's function, G is the grav- 
itational constant, and h(x) is an arbitrary vector 
field. Using the Green's function, equation J3J) can 
be written in the form 



M( — )a = gN 

ao 



V x h 



(5) 



where the first term on the right-hand side cor- 
responds to the Newtonian gravitational accelera- 
tion and results from the distribution of the matter 
density. Since, with respect to the right-hand side 
of equation ©, V ■ (V x h) — 0, we can conclude 
from Gauss's theorem that on a boundary at infin- 
ity, Vx/i falls faster than 1/r 2 , in other words, it 
can be ignore d in comparison with the Newt onian 
acceleration ( Bekenstein and Milgrom 1984h . For 
the interpolating function fJ,(x) given in equation 
@, we can obtain a in terms of gN as follows: 



a = 



l + (^£)2 

a N 



1/2 



(6) 



Here we model the Galactic disk as a simple, in- 
finitesimally thin Kuzmin disk (Kuzmin 1956), 
for which the corresponding Newtonian potential 
(Binney & Tremaine 1987) is given by 



4> N (R,z) 



—GM 



(7) 



where a is the disk length scale and M is the mass 
of the disk. The question can be asked whether we 
can express the acceleration in the MOND regime 
in terms of the gradient of a scalar potential. To 
answer this question, we write the MONDian ac- 
celeration in terms of the Newtonian potential 
from equation ((5]) as follows: 



a = V^Ar/(|V^Ar|). 



(8) 



The acceleration is compatible with a scalar poten- 
tial if the curl of the right-hand side of equation © 
vanishes. The condition V x [Vyjv/dVyjjvD] = 
implies that V/(|Vyjv|) x Wip^ = 0, so in the case 



3 



of a Kuzmin disk, for which |Vyjv| = l f%/GM, 
we have V/(|Vyj\r|) cx V(fN, which means that 
V x a — 0. A technique for constructing more 
general density-potential pairs in MOND has been 
introduced recently by Ciotti et al. (2006). 

2.2. CDM Halo Models 

In this subsection, we review the various halo 
models for the dark matter structure of Galaxy. 
Generic classes of these models are axisymmet- 
ric "power-law" and "logarithmic" models for the 
dark halo of the MW (Evans 1993, 1994). The 
gravitational potentials of the power-law and log- 
arithmic models are given by 



V 2 R ( , 



[3(R C + R 2 + z 2 q~ 2 )P/ 2 



(9) 



and 



<S>G = --V a 2 \og{R c 



R 2 + z 2 q- 2 ) (10) 



respectively. Poisson's equation gives the density 
distribution of the halo in cylindrical coordinates 
as 



p{R,z) 



V 2 R C 



(11) 



AirGq 2 

R 2 {1 + 2q 2 ) + R 2 {1- (3q 2 ) + z 2 [2 - (1 + f3)/q 2 
(R 2 + R 2 + z 2 /q 2 )^+ i )/ 2 

where the case (3 = corresponds to a logarith- 
mic model R c is the core radius, and q is the 
flattening parameter (i.e., the axial ratio of the 
concentric equipotential). The parameter deter- 
mines whether the rotation curve asymptotically 
rises (j3 < 0), asymptotically falls ((3 > 0), or is flat 
(/3 = 0). At asymptotically large distances from 
the center of the Galaxy in the equatorial plane, 
the rotation velocity is given by V C i rc ~ i? _/3 . The 
parameters of our halo models arc given in Table 
1. 

For the Galactic disk, we use for the density of 
the disk the double exponential function 



p{R,z) 



So 
2h 



exp 



R — i?o 






Rd 


exp 


h 



(12) 



(jBinnev fe Tremaine 19871 ). where z and R are 
cylindrical coordinates, Rq is the distance of the 
Sun from the center of the MW, Rd is the length 



scale, h is the scale height, and So is the column 
density of the disk. The parameters of the disk 
models are also indicated in Table 1. The Galac- 
tic bulge is considered to be a point-like structur e 
with mass (1 - 3) x 10 10 M Q (|Zhao fc Mao 19961) . 



3. The Magellanic Stream 

The MS is a narrow band of neutral hydro- 
gen clouds with a width of about 8° that extend 
along a curved path, starting from the Magel- 
lanic Clouds (MCs) and oriented toward the south 
Galactic pole (Van Kuilenburg 1972; Wannier & 
Wrixon 1972; Mathewson et al. 1974). This 
structure comprises six discrete clouds, labeled 
MS I to MS VI, lying along a great circle from 
(I = 91°, b = -40°) to (/ = 299°, b = -70°). 
The radial velocity and column density of this 
structure have been measured by many observa- 
tional groups. Observations show that the ra- 
dial velocity of the MS with respect to the Galac- 
tic center changes monotonically from Okms^ 1 
at MS I to -200fcms" 1 at MS VI, where MS 
/ is located nearest to the center of the MCs 
and the MW, with a distance of ~ 48 kpc, and 
MS VI is at the farthest distance. Another fea- 
ture of the MS is that th e mean column density 
falls from MS I to MS IV dMathewson et al-1977t 
iMathewson et al. 198TI ; iBriins et al. 20051 ). 



The MS is supposed to be an outcome of the 
interaction between the MCs and the MW, where 
the gravitational force of the MW governs the 
dynamics of the MS dMathewson et al. 1987 



Westerlund 1990HFuiimoto fc Murai 1984tlWavte 19891 ). 
There are several models for the origin and 
existence of the MS, such as tidal stripping 



dLin fc Lvnden-Bell 19821: iGardiner fc Nogiichi 1991 
IWeinberg 200d) . ram pressure (jMoore 1994HSofue 1994) 
and continued ram pressure stripping ( Liu 1992f h 
In the diffuse ram pressure model, there is a dif- 
fuse halo around the Galaxy that produces a drag 
on the gas between the MCs and causes weakly 
bound material to escape from the region to form a 
tail. The discrete ram pressure models are based 
on collision and mixing of the density enhance- 
ments of the halo with the gas between the MCs 
to form the MS, and in the scena rio of continued 
ram pressure stripping (Liu 19921) . the MS is in a 
state of quasi-steady flow (i.e., §f = 0). In the 
tidal model, the condition for quasi-steady flow 
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Model : 


S 


A 


B 


C 


D 


E 


F 


G 




Medium 


Medium 


Large 


Small 


Elliptical 


Maximal 


Thick 


Thick 


Description 


disk 


disk 


halo 


halo 


halo 


disk 


disk 


disk 


P 







-0.2 


0.2 














q 




1 


1 


1 


0.71 


1 


1 


1 


v a (kms~ l ) 




200 


200 


180 


200 


90 


150 


180 


R c (kpc) 


5 


5 


5 


5 


5 


20 


25 


20 


R (kpc) 
E o (Af pc- 2 ) 


8.5 


8.5 


8.5 


8.5 


8.5 


7 


7.9 


7.9 


50 


50 


50 


50 


50 


80 


40 


40 


R d (kpc) 


3.5 


3.5 


3.5 


3.5 


3.5 


3.5 


3 


3 


Zd(kpc) 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 


1 


1 



Table 1: Parameters of the Galactic models. The first line indicates the description of the models; the second 
line, the slope of the rotation curve ((3 — flat, j3 < rising and (3 > falling); the third line, the halo 
flattening (q = 1 represent spherical); the fourth line, (v a ), the normalization velocity; the fifth line, R Cl the 
halo core; the sixth line, the distance of the Sun from the center of the Galaxy; the seventh line, the local 
column density of the disk (E = 50 for canonical disk, E = 80 for a maximal thin disk and E = 40 for a 
thick disk); the eighth line, the disk length scale; and the ninth line, the disk scale hight. 



may be fulfilled by means of the tidal force of the 
Galaxy, which strips the HI clouds from the MCs. 
The Stream material then follows the same orbit 
as the MCs, such that the local velocity of the MS 
remains constant. 

Recently, Brims et al. (2005) reported on high- 
resolution observations of the radial velocity and 
column density of the MS. In the next section, we 
use these MS data to compare with the expected 
dynamics from the MOND and CDM scenarios. 
Comparing with the observations will enable us to 
place constraints on the parameters of the model. 

4. Dynamics of The Magellanic Clouds 

In this section, we obtain the dynamics of the 
MCs in MONDian theory, considering for the 
structure of the Galaxy an infinitesimally thin 
Kuzmin disk without a dark halo component. In 
this model, since we have ignored the Galactic 
halo, there is no drag force from the halo on the 
MS, and the main factor in the stripping of the 
MS is the tidal force exerted by the Galactic disk. 
Using the gradient of the gravitational force im- 
posed by the Galaxy on the MCs, we can make 
an estimate of the tidal radius, the distance from 
the MCs beyond which all the material can escape 
from the structure. 

In the Newtonian case, we consider a struc- 
ture with mass M s rotating in an orbit around 



the Galaxy, which has a halo of size Rg and 
mass Mq. The acceleration at the center of this 
structure is GM G / 'R g , while the acceleration at 
a distance r from its center toward the Galaxy 
is GM G /(R G - r) 2 , for R G » r, the differ- 
ence between these will be gc = 2GM G r / R G . 
On the other hand, the acceleration imposed to- 
ward the structure is g s — GM s /r 2 . So, mate- 
rial will move away in the direction of the Galaxy 
if g G > gs, which i mplies a stability radius of 
rudai = i? G (2^) 1 / 3 (jvon Hoerner 1957t >. For the 
deep MOND regime, fx{x) ~ x and the MONDian 
and Newtonian accelerations are related through 
a ~ (aoffjv) 1 ^ 2 - The tidal acceleration, similarly 
to the Newtonian case, obtains as the difference 
between the acceleration in the MCs and at a dis- 
tance r from them as 

g G = GM G aa 1/2 r/R G . On the other hand, 
we have the gravitational acceleration toward the 
MCs at the distance r, which is g s = GM s ao^ 2 /r, 
and so for the tidal radius we obtain T^dai — 
Rgijj^:) 1 ^ 4 - If we take the same mass-to-light ra- 
tio for the Galaxy as for the Clouds, we can say 
that the tidal radius in the MONDian regime is of 
the same order as that in the dark matter model. 
Let the ratio of the luminous matter of the MCs 
to that in the Galactic disk be of order 10 -2 ; for 
R G ~ 50kpc, we obtain in the MONDian regime 
a tidal radius of about ~ 15kpc, versus ~ Yikpc in 
the dark matter model. 
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Model 


s 


A 


B 


C 


D 


E 


F 


G 


MOND 


MOND (M. L.) 


x 2 


526 


816 


2285 


437 


576 


1511 


562 


315 


810 


356 



Table 2: Minimum x 2 from fitting with various Galactic models. The Cold Dark matter halo models are 
labeled with of <S, A, ... G with the corresponding values for the parameters of the models, indicated in 
Table [U The last column corresponds to the minimum % 2 from the Maximum Likelihood analysis in the 
MOND model. 



Model 


M(M ) 


a 


b 


c 


vo 


R c 


q 


M b (M ) 


ao(ms— 2) 


x 2 


MOND 


1.2 x 10 11 


4.5 














1.2 x 10" 1U 


356 


CDM 


1.2X10 11 


4.5 






175 


13 


l 






343 


J095CDM 


1.2 x 10 11 


4.5 






175 


13 


0.95 






349 


f09CDM 


1.2 x 10 11 


4.5 






175 


13 


0.9 






361 


fl.25CDM 


1.2 x 10 11 


4.5 






175 


13 


1.25 






341 


£05 


1 x 10 11 


6.5 


0.26 


0.7 


171 


13 


0.9 


3.53 x 10 10 




362 



Table 3: Comparison of MS data with the MOND and CDM Galactic halo models. The first line is the usual 
MOND model, the second to the fifth lines correspond to the CDM models (Johnston et al. 2005 and Law 
et al. 2005 ) The parameters a, b, c and R c are in kpc and vo is in terms of km/ 's. See Read and Ben Moore 
(2005) for more details of the models. 

drically symmetric distribution of matter around 
the z-axis. To solve the equation of motion, we 
take the location and velocity of MS / as the 
initial conditions for the dynamics of the MCs. 
Since MS I is located closest to the center of 
Galaxy, it almost has only a transverse veloc- 
ity component with respect to the Galactic cen- 
ter. Although, there are no indicators at MS / 
by which to measure its transverse velocity, we 
can nevertheless use the transverse velocity of 
the MCs measured from the kinematics of plan- 
etary nebulae in the LMC , vg = 281 ± 41fcms~ 1 
( van der Marel et al. 2002) as an estimate of the 
transverse velocity of MS /. We utilize the con- 
servation of angular momentum for the dynamics 
of the MCs (i.e., r M CsVMCs = tmsivmsi), which 
provides an MS / transverse velocity on the order 
of vg(ro) = 320 ± 50fcms _1 .Here the distance of 
MS / has been taken to be tq = 48 ± lkpc. 

In what follows, we compare the dynamics of the 
MCs with observation. Figure Q] shows the ob- 
served radial velocity of the MS in terms of the 
angular separation of from the MCs in a frame 
centered on the Galaxy. Since we are interested 
in comparing the global motion of the MS with 
the model, the internal dynamics of the structure 
has been ignored. The radial velocity of the MS 



According to the quasi-steady flow model, the 
MS will lie in the same orbital path as the MCs, so 
the location and velocity of the MS can represent 
the dynamics of the MCs. Here we are looking to 
the radial velocity of the MS toward the Galactic 
center to compare with observations. As a result of 
hydrodynamic friction, the MCs could have their 
motion damped and fall toward the Galaxy in a 
spiral path. However, since in this model there is 
no halo around the Galaxy, we can ignore friction 
forces. Substituting the Newtonian potential of 
the Kuzmin disk (eq. [7]), the MONDian potential 
in cylindrical coordinates is obtained 



(j)M= VGM^ HR2 + {a+Mn (13) 
as (|Read fc Moore 20051: iBrada fc Milgrom 19951 ). 



For convenience, we carry out a coordinate trans- 
formation from the cylindrical to the spherical 
frame and obtain the components of the accelera- 
tion as follows: 

a cos 8 + r 



■y/MGaoi- „ 2 14 

2ar cos + r l + a 1 



ag 



2ar cos + r l + a z 
fl0 = 0, (16) 

where is zero as a consequence of the cylin- 
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Angular distance (radian) 



Fig. 1. — Radial velocity v r vs. angular dis- 
tance 9 along the MS in the center of Galacto- 
centric frame. The observational data (Brians et 
al. 2005) are compared with theoretical curves 
obtained from MONDian and CDM models. The 
minimum x 2 correspondsing to each curve is given 
in TablelSJ^ee the electronic edition of the Journal 
for a color version of this figure.] 



is calculated by averaging over the data at each 
point of the Stream and the error bars on the ra- 
dial velocity result from the dispersion velocity of 
the structure. This velocity dispersion is caused 
by substructural motions of the MS such as the 
stochastic behavior of gas, and for different real- 
izations of the MS, the dispersion velocity will be 
different at each point. The observed data are a 
continuum distribution of radial velocities in terms 
of angular separation; however, we collected data 
for each cluster of gas and obtained the average 
velocity and the corresponding dispersion veloc- 
ity. We use x 2 fitting, 

N 

x 2 = E( ) 2 > 

l 

for comparison of the observed data and dynam- 
ics of the MS from the MOND and CDM models 
(see Fig. [TJ). We should take into account that 
the CDM and MOND models in these cases have 
different mass distributions. 




a 

Fig. 2. — Marginalized likelihood functions of the 
parameters of the MOND and Kuzmin disk. Top, 
product of the disk mass Mq (10 11 Mq) and the 
MOND acceleration scale ao (10~ 10 ms~ 2 ); bot- 
tom, length scale of the Kuzmin disk. The in- 
tersection of the curves with the horizontal solid 
and dashed lines give the bounds at la and 2a 
confidence levels, respectively. [See the electronic 
edition of the Journal for a color version of this 
figure.) 
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id 




M „ a. 

Fig. 3. — Joint confidence intervals for 
M o a o (lO 1:L M0 x lCr 10 TOs- 2 ) and the disk length 
scale a(kpc), with la (solid line), 2a (dashed 
line), and 3<r long-dashed line confidence levels. 
The minimum value of \ 2 corresponds to a = 
0.2j£f kpc and Af a = 0.84±g;g|. [See the elec- 
tronic edition of the Journal for a color version of 
this figure.] 
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Fig. 4. — Constraints on Moao = 8.4+q '\ horn the 
maximum likelihood analysis, with la confidence 
interval (M in terms of 10 10 M Q and MONDian 
acceleration scale ao in terms of 10~ 10 ms~ 2 ). Us- 
ing the range ao = 0.8 — 1.4 implies Mo = 6 — 11 
for the mass of the Kuzmin disk. [.See the electronic 
edition of the Journal for a color version of this 
figure.} 



The dark matter models are labeled S, A, B, C, 
D, E, F, and G with their corresponding parame- 
ters in Table [TJ The minimum x 2 values from fit- 
ting with the CDM and MOND models are listed 
in Table [21 We performed a maximum likelihood 
analysis to find the best parameters for MOND. 
Comparing the minimum value of \ 2 from the 
modified dynamics with those of the dark mat- 
ter models shows that MOND results in a bet- 
ter fit than the CDM models, except for model 
G. The best parameters for a MONDian poten- 
tial and the characteristic acceleration parameter, 
from the maximum likelihood analysis, are found 
to be M a = 0.84±g;g!j and a = 0.2±g;f 5 kpc, 
where Mq is the mass of the disk in terms of 
1O 11 M0, ao is the MOND acceleration scale in 
terms of 10~ 10 ms~ 2 and a is the length scale of 
the Kuzmin disk. The marginalized relative likeli- 
hood functions for MqOq and a are shown in Fig- 
ure[2j Figure [3] shows the joint confidence intervals 
for Moao & n( i a with la to 3a confidence levels. 



As can be seen from the results of the maximum 
likelihood analysis, in MOND the length scale of 
the Kuzmin disk favor small values, in the range 
a < lkpc with a wide error bar. We can interpret 
this as a consequence of the MS being located at a 
large distance from the Galaxy compared with a, 
which means that the dynamics of the MS is not 
sensitive to this length scale. Moao is also con- 
fined to a range of about 0.8 to 0.9. In order to 
break the degeneracy between the mass of the disk 
and the acceleration scale of MOND (ao), we use 
the luminous mass of the MW to constrain the pa- 
rameter ao- The total luminosity of the Galacti c 



disk is 1.2 x 10 10 L Q (jBinnev fc Tremaine 1987h . 



and if we assume an average stellar mass-to-light 
ratio of T = 5, the total stellar mass of the Galac- 
tic disk should be about M* ~ 6 x 10 10 M Q ; using 
the constraint Afoao — 0.84 then results in ao — 
1.4. Comparing the value ao = 1.2 ± 0.27 from 
analysis of the rotati on curves of spiral galaxies 
(jBegeman et al. 19911 ). with that from our analy- 
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Fig. 5. — Constraints on Mqcio considering 2a er- 
ror bars on the initial value for the transverse ve- 
locity of MS /. Using the preferred value M* ~ 6 x 
1O 1O M for the Galactic disk yields a = 0.8-1.4; 
alternatively, adopting a MOND acceleration scale 
in the range ao = 0.8 — 1.4 implies a disk mass in 
the range M = (2.5-23) x 10 10 Af Q .[5"ee the elec- 
tronic edition of the Journal for a color version of 
this figure.] 



sis of the MS, we see good agreement. If, in addi- 
tion to the contribution of the disk to the Galactic 
mass, we add the lower and upper limits for the 
mass of the Galactic bulge, the total mass of the 
Galaxy lies in the range (4.3 - 12.8) x 10 10 M m 



([Sackett 1997t IZhao et al. 19951) . which confines 



the value of ao to 0.65 — 1.95. In addition, for the 
range of the MOND acceleration scale adopted in 
the literature, ao = 0.8 — 1.4, the mass of the disk 
is found to be M = (6 - 11) x 10 10 M Q (see Fig. 

a. 

In the likelihood analysis here, we used as an 
initial condition the central value of the velocity 
of MS I without taking into account the corre- 
sponding error bar. In order to check the sen- 
sitivity of our results to the uncertainty in the 
velocity of the MCs, we repeated the likelihood 
analysis for an MS I transverse velocity in the 
range vmsi — 320 ± 100 with 2a error bars. Fig- 
ure [5] shows that Mqclo is constrained to MqOq = 



Fig. 6. — Same as Fig. 1, but for recent CDM 
models of the Galactic halo form studies of the 
tidal debris from the Sagittarius dwarf galaxy (see 
Table [3|). [See the electronic edition of the Journal 
for a color version of this figure.] 



0.84^^53 when considering this uncertainty on the 
initial value for MS I's transverse velocity. Us- 
ing the preferred value of M* ~ 6 x 10 10 Mq for 
the Galactic disk confines the acceleration scale to 
the range a = 0.8 — 1.4. Alternatively, consider- 
ing the literature acceleration scale for MOND of 
ao = 0.8 = 1.4 provides a Galactic disk mass of 
M = 2.5 - 23 x 10 10 M Q , a larger range than that 
of our previous analysis. For the standard value 
of ao = 1.2, the mass of the Galactic disk is con- 
fined to Mq = 3 — 15 x 10 10 Af Q . Higher accuracy 
observations of the transverse velocity of the MCs 
will provide a better constraint on the parameters 
of the model. 

Finally, we compare the dynamics of the MS 
in MOND with recent CDM Galactic halo models 
proposed to interpret the dynamics of the Sagit- 
tarius dwarf galaxy around the MW (Read & 
Moore 2005). Figure [HI compares the observed ra- 
dial velocity of the MS with these Galactic models, 
and Table [3] indicates the best-fit value of x 2 f° r 
each model and shows that the MOND result is 
almost compatible with the CDM models. 
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5. Conclusion 

In this work, we used recent observational data 
on the radial velocity of the Magellanic Stream 
to examine modified Newtonian dynamics. The 
MS is considered to be a trace of the Magel- 
lanic Clouds produced by the interaction of the 
Clouds with the Galaxy. We considered the 
MONDian tidal effect of the Galaxy as the ori- 
gin of the MS. For simplification, in our analy- 
sis we used the quasi-steady model for the flow 
of the MS, considering this structure to follow 
the same orbit, with the same dynamics, as the 
MCs. The radial velocity of this structure was 
compared with the observations, allowing us to 
place constraints on the product of the Galac- 
tic disk mass and the MOND acceleration scale 
(Motto). The range for the mass of the Galaxy 
obtained from this analysis is compatible with 
the observed luminous mass of the disk. Com- 
paring the best fit from MOND with cold dark 
matter models for the halo shows that MOND 
fits the obser vations as well as the conventional 
halo models dRead fc Moore 2005t iHelmi 2004 



Johnston et al. 2005 ; Hasani-Zonooz et al. 20051 ). 

We should point out that while the dynamics 
of the MS in this model is in good agreement 
with the observations, a problem with this tidal 
model could be the lack of corresponding stel- 
lar tidal debris in the MS fdMoore fc Davis 1994 : 



Guhathakurta & Reitzel 19981) ). Despite our ex- 
pectation in the tidal model that stars should have 
been stripped off with the gas, observations show 
that stellar tails are often completely missing or 
are offset from the gaseous tails in interacting sys- 
tems (e.g., the M82 Group; Yun et al. 1994; Hi- 
bbard et al. 2000). In the case of a MONDian 
tidal model, since the tidal radius is larger than in 
the Newtonian case, we expect a Galactic struc- 
ture with a gas distribution that is more extended 
than the stellar component form of the MS. So, we 
should be observing mainly the gaseous compo- 
nent in the MS. N-body simulations of tidal strip- 
ping of the MS from the MCs in the MONDian 
regime without a Galactic halo will give a better 
view of this model and enable us to compare the 
density distribution of the gas in the MS with ob- 
servations. 



manuscript and providing very useful comments, 
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Institut d'astrophysique de Paris, where a part of 
this work was done, for their hospitality. 

REFERENCES 

Alcock, C, et al. 2000, ApJ, 542, 281 

Baumgardt H., Grebel E.K. and Kroupa P. 2005, 
MNRAS, 359, LI 

Bekenstein J.D., Milgrom M., 1984, ApJ, 286, 7 

Bekenstein J., 2004, PRD, 70, 083509 

Begeman K.G., et al, 1991, MNRAS, 249, 523 

Binney S., Tremaine S., 1987, Galactic Dynamics, 
Princeton Univ. Press, Princeton, NJ 

Brada R., Milgrom M., 1995, MNRAS, 276, 453 

Brims, C, et al. 2005, A&A, 432, 45 

Buries, S., Nollett, K. M., & Turner, M. S. 2001, 
ApJ, 552, LI 

Ciotti L., Londrillo P. and Nipoti C, 2006, ApJ, 
640, 741. 

Evans N. W., 1993, MNRAS, 260, 191 

Evans N. W., 1994, MNRAS, 267, 333 

Famaey B., Binney J., 2005, MNRAS, 361 633 

Fujimoto M., Murai T., 1984, Structure and Evo- 
lution of the Magellanic clouds, IAU symp. No. 
108, 115. 

Gardiner, L. T., Noguchi, M. 1996, MNRAS, 278, 
191 

Goodman M., Witten E., 1985 , Phys. Rev. D31, 
3059. 

Green A.M., 2003, to appear in the proceedings of 
IAU Symposium 220 "Dark matter in galaxies" , 
ASP, Eds: S. Ryder, D.J. Pisa no, M. Walker, 
K. Freeman, (|astro-ph/0310215 ) 



The authors thank J. I. Read for studying the 



10 



Guhathakurta, P., Reitzel, D. B. 1998, in ASP 
Conf. Ser. 136, Galactic Halos: A UC Santa 
Cruz Workshop, ed. D. Zaritsky (San Francisco: 
ASP), 22 

Hasani-Zonooz A., Haghi H., Rahvar S., 
flastro-ph/050417lD 

Helmi A., 2004, ApJ, 610, L97 

Hibbard, J. E., Vacca, W. D., Yun, M. S. 2000, 
AJ, 119, 1130 

Johnston K.V., et al , 2005, ApJ, 619, 800 

Jungman G., et al, 1996, Phys. Rep. 267(no 5-6), 
195 

Kuilenburg , J. V., 1972, A&A, 16, 276 

Kuzmin , G. G., 1956, AZh, 33, 27 

Kirkman, D., Tytler, D., Suzuki, N., OMeara, J., 
& Lubin, D., 2003, ApJS, 149, 1 

Lin D. N. C, Lynden-Bell D., 1982, MNRAS 198, 
707. 

Liu Y., 1992, A&A 257, 505 

MathewsonD.S., ClearyM.N., Murray J.D., 1974, 
ApJ, 190, 291 

Mathewson D.S., Schwarz M.P., Murray J.D., 
1977, ApJ, 217, L5 

Mathewson D.S., Wayte S.R., Ford V.L., Ruan K., 
1987, Proc. Astron. Soc. Aust. 7, 19. 

Meatheringham S. J., et al, 1988, ApJ, 327, 651 

Milgrom M., 1983, ApJ, 270, 365 

Mirabolfathi N., Invited talk at the XXIV Physics 
in Collisions Conference (PIC04), Boston, USA, 
June 2004, ( |astro-ph/0412T03| ) 

Moore B., 1994, Nature 370, 629. 

Moore B., Davis M., 1994, MNRAS 270, 209 

Nusser, A., Pointecouteau, E., 2006, MNRAS, 
366, 969 

Paczyhski B., 1986, ApJ, 304, 1 

Pointecoute au E., Silk P. , 20 05, MNRAS, 364, 
654-658, ( |astro-ph/05050T7| . 



Pretzl K., 2002, Space Science Reviews, v. 100, 
Issue 1/4, p. 209-220. 

Rahvar S., 2004, MNRAS, 347, 213. 

Read J.I , Ben Moore, 2005, MNRAS, 361 , 971- 
976. 

Sahu K, C, 2003, "To appear in the Pro- 
ceedings of the STScI Symposium on "Dark 
Universe: Matter, Energy, and Gravity" 
( |astro-ph/0302325[ ). 

Sanders, R. H., 2001, ApJ, 560, 1 

Sanders R. H., McGaugh S., 2002, Ann. Rev. As- 
tron. Astrophys., 40, 263. 

Sackett P.D., 1997, ApJ, 483, 103-110 

Scarpa R., Marconi G., Gilmozzi R., 2003, 
IAU220 symposium "Dark matter in galaxies", 
( |astro-ph/0308416[ ) 

Sofue, Y., 1994, PASJ, 46, 431 

Spergel, D. N., et al, 2003, APJS 148, 175 

Tisserand, P., Milsztajn, A., 2004, appear in the 
proceedings of the 5th Rencontres du Vietnam, 
" New Views on the Universe" , 

van der Marel R.P., et al, 2002, AJ, 124, 2639-2663 

von Hoerner S., 1957, ApJ, 125, 451 

Wannier P., Wrixon G.T., 1972, ApJ, 173, L119 

Wayte S. R., 1989, Proc. Astron. Soc. Aust.8, 195 

Westerhmd B. E., 1990, A&A R 2, 29 

Weinberg, M, D., 2000, ApJ, 532, 922 

Yun M. S., Ho P. T. P., Lo K. Y. 1994, Nature, 
372, 530 

Zhao H., Mao S., 1996, MNRAS, 283, 1197Z. 

Zhao H., Spergel D.N., Rich R.M., 1995, 
ApJ,440L, 13Z 

ZhaoH., Bacon D. J., T aylor A.N., Horne K., 2006, 
MNRAS, 368, 171 ( astro-ph/0509590| ). 

Zh ao H.S, Famaey B, 2006, ApJ,638, L9-L12 
( |astro-ph/0512425D 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 



11 



